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The excitation by a high frequency field of multi-level 
quantum systems with a slowly varying density of states is 
investigated. A general approach to study such systems is 
presented. The Floquet eigenstates are characterized on sev- 
eral energy scales. On a small scale, sharp universal quasi- 
resonances are found, whose shape is independent of the field 
parameters and the details of the system. On a larger scale 
an effective tight-binding equation is constructed for the am- 
plitudes of these quasi-resonances. This equation is non- 
universal; two classes of examples are discussed in detail. 

The excitation of quantum systems by external fields, 
is a fundamental problem in quantum mechanics. Any 
interaction between matter and electromagnetic fields is 
described, on the simplest level, by such a model. Given 
the system at some initial condition, one would like to 
characterize its energy absorption from the external field 
as a function of time. In some cases, simple quantum sys- 
tems can be described effectively by one degree of free- 
dom. One-dimensional (ID) driven systems have been 
studied widely within the field of "quantum chaos" , 
since they provide some of the simplest examples of quan- 
tum systems which are chaotic in the classical limit. The 
phonomenon of dynamical localization, where the classi- 
cal energy absorption associated with chaotic motion is 
suppressed by quantum interference effects || , has been 
of special interest (For reveiew see S. Fishman, in (I)). 

In this work, a general approach for systems with 
a slowly-varying density of states, driven by a high- 
frequency smooth periodic field, is presented. We study 
the quasi-energy (Floquet) eigenstates, which are the 
stationary states of the time dependent system B, in 
the unperturbed representation, and characterize their 
properties over several energy scales. On the smallest 
scale, we find that, under quite general conditions, these 
states are composed of a ladder of sharp peaks, or "quasi- 
resonances" (QR). These are related to quantum nearly 
resonant transitions between the energies of the undriven 
system and should be distinguished from the classical res- 
onances. While the existence of such peaks is well known 
||, here they are derived from a local exact solution, 
which enables the quantitative description of their loca- 
tion and shape. Surprisingly, the QR's turn out to have 
a universal shape, independent of field parameters and 
of the details of the system. In particular, their width is 
independent of the driving field strength, a result which 
contrasts simple perturbation theory or a two-level ap- 
proximation. On a larger energy scale, our local solution 
allows us to construct an effective equation for the enve- 
lope of amplitudes superimposed on the QR's ||. The 



form of the local solution is used explicitly to obtain a 
tight-binding equation on the lattice of QR peaks. The 
parameters of this equation are non-universal and are re- 
lated to the spectrum of the unperturbed system and to 
the matrix elements of the perturbation. 

Investigation of the Floquet operator (the evolution 
operator for one period), is very convenient for time pe- 
riodic systems. Its spectrum determines the quantum 
dynamics on all time scales, similar to the Hamiltonian 
for autonomous systems. Although finding this spec- 
trum exactly is as difficult as solving the time-dependent 
Schroedinger equation, here we characterize the eigen- 
states by general qualitative properties, from which prop- 
erties of the dynamics follow. This idea is inspired by the 
kicked rotor analogy R, where a robust qualitative fea- 
ture of the quasi-energy states - namely their exponen- 
tial localization in energy space - explains the quantum 
saturation in energy absorption. A formal mapping was 
constructed to a ID Anderson model, which relies on spe- 
cial properties of the kicked rotor; other driven systems 
require in general a different treatment. 

For systems with a slowly varying density of states, 
a constant energy spacing is well defined on a local en- 
ergy scale, and the deviations of the spectrum from har- 
monicity appear as an adiabatic change of this spacing 
with energy. This "adiabatic nonlinearity" of the spec- 
trum is characteristic of many ID systems, of which sev- 
eral examples are the hydrogen-like atom ( (lj], D. L. 
Shyepelvansky in [[lj), charge bubbles in liquid helium 
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and surface electrons in a 2D metal in a per- 
pendicular magnetic field This work presents a gen- 
eral framework for this class of systems, from which some 
new results emerge and some known results for special 
cases are confirmed. We focus on a monochromatic driv- 
ing, which is of high-frequency compared to the typical 
frequencies of the system. Thus, many unperturbed lev- 
els participate in the excitation (it is important, how- 
ever, that the spectrum is discrete for our solution to 
hold). The basic idea employed in this work, is to solve 
the problem in a limited energy regime, where the locally 
defined energy spacing can be considered constant. This 
local solution relies on the exact quantum mechanical so- 



lution of an integrable model |13|. We then account for 
the large energy scales by exploiting the adiabatic de- 
pendence of parameters on energy. A similar idea was 
previously applied to the "bubble" model [Q. Here a 
wider class of systems is studied (more details in [^s|). 

Consider the following ID Hamiltonian in action-angle 
variables (I, 9) of the bound system: 



H 



H (I) + k V(I)g(0)cos(Slt). 



(1) 



It is assumed that TLo(I) and V(I) are smooth func- 
tions of /, and that g{9) has a Fourier expansion with 
smooth, slowly decaying components G(m). Our object 
is to calculate the Floquet operator, and to character- 
ize its eigenstates in the unperturbed representation |n), 
where / = nh. In the high-frequency regime, the field 
period is a relatively short time scale and thus the semi- 
classical approximation is expected to be very accurate. 
Using also a leading order approximation to the classical 
trajectories |l(|, the Floquet eigenvalue equation is 
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where T = 2n/Q is the period of the external driving, A 
is the quasi-energy, and 

A(I,6)=k [ V(I)g(9 + u(I)t)cos(ni)dt. (3) 
Jo 

(for details, see The physical reason why the Flo- 

quet operator can be factored into a product of two op- 
erators, as in Eq. (Q) , is that most of the energy transfer 
takes place around the singular point of the potential. 
When approximating a 3D potential by one degree of 
freedom, usually one is dealing with the half line so a 
special point naturally arises. For the hydrogen atom, 
e.g., this is the nucleus; this has been the basis for con- 
structing the Kepler map |l7]|n| • The existence of a sin- 
gular point in space is related to the slow decrease of the 
Fourier components of the driving, G(m). 

We define a dimensionless parameter e = uj/Q, where 
u> = dHo/dl. In the absence of the driving, the solutions 
of Eq. (||) are 5-functions in n-space. For e <C 1, the per- 
turbation couples most effectively (^-functions separated 
by approximately jhVL in energy (j integer), correspond- 
ing to nearly degenerate A's. These form the QR ladder, 
which will be derived in what follows. 

For e < 1, 

A(I;6) = { 0<e<2«(l-e) 
V ; \ cos[^] 2^(l-e)<0<27r V ' 

where R — =£-kV{I)G{m+), and m + is the integer closest 
to 1/e, corresponding to classical resonance = m + u. 
The Fourier series of A is dominated by this resonance; 
the smoothness of the driving coefficients allows us to 
take into account not only this term but also its vicinity. 

Since 7io(nh) varies slowly as a function of n, it can 
be expanded around a large value hq to first order: 
Ho{nfi) ~ Hq{jioK) +u>IK, where I = (rt — no). Then the 
eigenvalue equation (^) appears as the equation for a "lin- 
ear kicked rotor" [|l3| , which is an exactly solvable model. 
Denoting the quasi-energies of the linearized model by 
X hn , its spectrum is \ hn T = [iu;T(mod 2n), and the cor- 



responding eigenstates are \u)e 
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, where 



J m (B) sinc[7r(^ — I + m/e)] (5) 



with the following definitions: 

<p = n(ji- 0(1 - 2e) 

B = £-kV(n h) G(m+)/sin(7r/e) 



(6) 



and sinc(:r) = sm(x)Jx. In the high-frequency limit, 
e<1, the function (g) is composed of a chain of peaks 
separated by approximately the energy of one photon, 
each weighted by an amplitude. These quasi-resonances 
(QR) are described by the function 
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sin[7r(Z — lj+Sj)] 
ir(l — lj+5j) 



(7) 



where lj is the center of the QR and 5j characterizes its 
precise shape. The Q/s constitute a nonorthogonal set 
which is probably complete, however for a given quasi- 
energy, the corresponding eigenstate (||) is composed only 
of a small fraction of the Qj's, at positions determined 
by \ Un through the relation E Xj = X^h+jfifl-TkjSj, with 
\6j \ < 1/2. Assuming that at each point tails of the sine 
functions centered far-away contribute incoherently, the 
absolute value of the wave function is large on this ladder 
of QR's. The width of the QR is independent of the 
driving field strength and of the density of unperturbed 
states. This is in contrast to the width associated with 
the transition rate given by Fermi's golden rule, and also 
in contrast to the Rabi width, both of which are invalid 
approximations in our regime of parameters. 

The next step is to construct the eigenstates of the 
original nonlinear system by matching different energy 
regimes, in each of which a local solution holds. In this 
construction, there are two important scales in action 
space. The first is the distance between classical reso- 
nances: our local solution is valid only outside the close 
vicinity of the classical resonances. The second is the 
width of the resulting local solution, which is determined 
by the Bessel function in Eq. (^) . Once the linearized lo- 
cal functions become wide enough to cover more than one 
classical resonance, it is exptected that the different local 
solution may interact and matching between the regions 
may be employed. It turns out that for the two classes 
of models described below, this matching condition co- 
incides with the Chirikov critetion for resonance overlap 
p9| . This implies that for these cases, the matching is 
valid inside classically chaotic connected regions of phase 
space. Outside such regions, the decay of the eigenstates 
is expected to be determined by classical bounds. 

Under conditions that the matching is valid, we use 
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(dE n /dn) no I, 



(8) 



to find that in the matched nonlinear solution, the QR 
appear on the ladder of states rij satisfying 



E,, 



\h + jhn-(dE n /dn) no 5 3 . 



(9) 



Thus the eigenvalue A sets the origin of the ladder for the 
corresponding eigenstate, and the QR's are numbered by 



their position j on this ladder. Each is described approx- 
imately by the function Qj (n) of Eq. (Q) , characterized 
by the peak position rij and the detuning \Sj\ < 1/2, both 
satisfying Eq. (^|). The matching may alter the am- 
plitudes of the QR's, therefore we write for the global 
eigenstate the following approximate form: 



(n\u) = ^2,A j Qj{n) 



S^ A . sin[7r(n-n J +^)] 
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nin-rij+Sj) 



This equation, together with (|^) which defines the pa- 
rameters rij and Sj, constitute the main result concern- 
ing the structure of the Floquet eigenstates on a small 
scale. They predict, for a given A, the precise location 
and shape of the QR's of the corrresponding eigenstate. 

We now turn to determine the amplitudes Aj. The 
Floquet eigenvalue equation can be equivalently written 
in an extended Hilbert space of variables (n, j) of which 
\n,j) = |n)e _l:,nt are basis functions (2(J: 



(E n - hjCl)</> n j + 



h\(j) n 



(11) 



where <j)n,j — (n,j\(f>) and O — V{I)g{9). For smooth 
functions V(I) the semiclassical matrix elements may be 
written as (n\(D\n') V(I)G(n — n'). Using the local 
structure of the QR's and the properties of the sine func- 
tions on a lattice, the equation reads: 



2hu 



nkV(hn,j) 



sin( 7 r^)A J - G + A j+1 - G~ Aj- X ~ (12) 



where G ± — G(j^- + Sj±i), and ej±i is the value of e at 
the (j ± l)-th QR. This is a tight-binding equation on 
the ID lattice of QR's labelled j. The diagonal potential 
is determined by the Sj, which are the detunings of the 
energies E n . from exact resonance, and the hopping is 
related to matrix elements of the driving between neigh- 
boring QR's. The fact that there is only near-neighbor 
coupling is a result of the driving being harmonic. 

We consider two classes of models in some detail. 
These are described by the Hamiltonian 
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k Ocos(fX), 



x > (13) 



where a is a real number, — 2 < a < 2. The driving field 
is described in the dipolc approximation, and it is conve- 
nient to work in different gauges for the cases of positive 
or negative a, thus O = x for a > 0, and O — p/fl for 
a < 0. Both classes have a slowly varying density of 
states, the interaction with the field is described by a 
function V(I) which is slowly varying with /, and g{6) 
with smooth Fourier coefficients. For high energies, ( |13| ) 
is well approximated by (|lj). Thus, on the fine scale they 
exhibit the universal structure of the QR's. The differ- 
ences between the models appear in the envelope of am- 
plitudes superimposed on the peak structure. 



For a > 0, the binding potentials have a triangle-like 
singularity near the origin x = 0. The special case of 
(7 = 1 corresponds to the bubble model, which is a tri- 
angular potential well. The dipole matrix elements have 
the asymptotic form (n\x\m) cx (nh) 2 ^ 2 ^ /(n — m) 2 . 
The eigenstates for this case may be characterized by two 
qualitatively different regimes in n-space. In the small n 
regime they are exponentially localized, while for larger 
values of n they are more extended, with a crossover point 
n c between the two regions, satisfying 



(n 2 /k) 



(2+<j)/<j 



(14) 



This is a generalization of a result previously found for 
the special case of the bubble model |l4|] . In the asymp- 
totic regime n — > oo, Eq. (|l2|) becomes similar to the 
Anderson model in a static electric field: 



sin(7rjj) \ ( Vnn 3 / 2 \ 
-7T) \-7hk-) Aj 



R + A jA 



FT A 



j'-i 



(15) 



where C\ is a constant and Br — (1 =F <5jiiej±i) 2 ■ The 
diagonal potential is given by the pseudo-random func- 
tion sm(nSj), where Sj = frac |i [ha(\ + j)} 1 ^^ and 

a = 2cr/(2 + cr), damped along the lattice by l/y^J. The 
hopping terms depend weakly on position j. For a ran- 
dom potential and strictly constant hopping, the eigen- 
states are power localized with a power proportional to 
the square of the prefactor of the diagonal potential [^l| . 
Numerical calculations for the model ([t5]) with (7 = 1, 
provide evidence that these differences do not alter the 
ualitative conclusion from the theorem. ( JI3], following 
) . Thus the eigenstates are power-law decaying, 
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(hn 3 /k 2 
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This im plies the existence of a crossover field strength 
k c ~ VM7 3 beyond which the tails of the eigenstates 
turn non-normalizablc. A similar result was stated pre- 
viously in H . For the special case of the bubble model, it 
was first obtained by use of a Kepler-like map |I1J . The 
reason that in this regime the result is independent of a 
is related to the dipole matrix elements being asymptoti- 
cally similar for all a, which is a result of the triangle-like 
singularity that all these potentials have at the origin. 

For a < 0, the potential is singular at the origin 
(linXj^o x a = — oo). The special case of a = — 1 corre- 
sponds to the ID Coulomb potential, which is a good ap- 
proximation for hydrogen-like rydberg atoms with small 
angular momentum ]22| . For this case, the derivation of 
the universal QR structure is strictly justified only as a 
leading order in the field strength k. Using the asymp- 
totic form (n\p\m) oc i sgn(n-m)(n?i)' T / afcr /|n — m| 2 / 2 ^, 
the resulting equation for the amplitudes is, 



sin(7r<5 r ) A r + R + A r+ i - i?~A r _i w (17) 



where C'2 is a constant, i? ± = (l=p<5 r±1 e,-±i) °- 2 , 77 = (4 — 

cr)/(2 - a) and <5 r = fracji [a (h\~ (r max -r)hQ)} 1/a ^ 

for r < r maa; ; since the bound part of the spectrum is 
of finite range in energy, there is only a finite number 
of QR's. For e < 1, the diagonal potential may be con- 
sidered similar to a random one. Neglecting the weak 
dependence of the hopping terms on r, one obtains a ID 
Anderson model for the amplitudes A r . The eigenstates 
of this model are exponentially decaying with a pertur- 
bative estimate for the localization length given by 



£ ~ k 2 /h 2 n 2r >. 



(18) 



Also in this case, numerical calculations confirm that the 
weak dependence of the hopping terms on position does 
not alter the existence of exponential localization and the 
scaling (|l8|). Our description neglects completely effects 
of the continuum, and thus is valid only far away from 
it; the result ( p"8| ) is meaningful only if £ is smaller than 
the lattice size, which is finite in this case. 

In conclusion, for a particle in a ID potential well, 
where the unperturbed spectrum is slowly varying, driven 
by a harmonic high-frequency field, some properties of 
the Floquet eigenstates were found on several energy 
scales. The small scale resonant structures, the QR's, are 
non-perturbative and universal, independent of the de- 
tails of the system and on the parameters of the driving. 
Thus the general form of the eigenstate is as described 
by Eqs. (0), ([To|). On a larger scale, it was shown in gen- 
eral how to construct a tight -binding equation, Eq. (12), 



for the amplitudes Aj of (|1C|) superimposed on the QR's. 
Two specific classes were considered, where this equation 
turns out to be an Anderson model on a finite lattice, and 
an Anderson model on an infinite lattice with a constant 
electric field. It will be of great interest to study exper- 
imentally the nature of the eigenstates and in particular 
the shape of the quasiresonances for some systems that 
belong to classes explored in the present work. Recent 
developments in the optical trapping of atoms may en- 
able such investigations for various potentials. 
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